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> ■ Abstract 

O . 

A rigourous theory for the determination of the van der Waals interactions in colloidal systems 
is presented. The method is based on fluctuational electrodynamics and a multiple-scattering 
method which provides the electromagnetic Green's tensor. In particular, expressions for the 

' 5-T, 

Green's tensor are presented for arbitrary, finite, collections of colloidal particles, for infinitely 
periodic or defected crystals as well as for finite slabs of crystals. The presented formalism allows 
for ab initio calculations of the vdW interactions is colloidal systems since it takes fully into account 
retardation, many-body, multipolar and near-fields effects. 

O . PACS numbers: 12.20.Ds, 42.50.Vk, 78.67.Bf 
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I. INTRODUCTION 



The van der Waals (vdW) interactions are particularly important in colloidal systems 
since, along with the electrostatic forces, they determine the structure of such systems. The 
stability of colloidal systems resulting from the interplay between the vdW and the elec- 
trostatic interactions is very well elucidated in the context of Derjaguin-Landau-Verwey- 
Overbeek theory- The vdW interactions which originate from the irreducible electro- 
magnetic (EM) fluctuations of vacuum are usually calculated by means of the Hamaker 
approach,- where the force stems from simple pairwise addition of the corresponding inter- 
molecular forces,- 1 ^ although the vdW interactions are not additive. A rigorous treatment 
of the vdW interactions based on ffuctuational electrodynamics*^ has been pioneered by 
Lifshitz^ for the case of two infinite half-spaces. The Lifshitz theory has been extended to the 
case of pairs of finite-sized objects such as spheres or cylinders (Derjaguin approximation)^ 
which is valid, however, for very short distances between the objects, in the nonretarded 
limit. In some cases, elements of the Lifshitz theory for half-spaces are incorporated within 
the Hamaker formula for the vdW force between two particles, in the form of semi-empirical 
corrections.-^&^ii By use of perturbation theory and the Clausius - Mossotti formula, 
Langbein^^ developed a general formalism for the vdW force between two spheres which 
has been primarily applied to aerosol particles.— dkd!L 

Recently, a new, rigorous theory based on ffuctuational electrodynamics for the calcu- 
lation of the vdW interactions among a collection of macroscopic bodies of finite size has 
been proposed.— This theory is based on a multiple-scattering Green's tensor formalism 
incorporated within the framework of fluctuational electrodynamics. More specifically, the 
vdW force results from the integration over the surface of the bodies of the Maxwell stress 
tensor of the vacuum/ thermal EM field which is provided by the fluctuation-dissipation 
theorem and through this by the Green's tensor of the classical EM field. The calculation of 
the Green's tensor is based on an EM multiple scattering formalism for arbitrary collections 
of scatterers. The multiple-scattering Green's tensor formalism offers a precise knowledge 
of the fluctuating EM field by going beyond the approximation of pairwise interactions 
between the scatterers and by taking into account the full multipole interactions between 
them. Furthermore, since it constitutes a solution to the inhomogeneous wave equation, 
retardation effects are included a priori in the presented formalism. In addition, metallic 
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and dielectric particles are treated on an equal footing since the method in question also 
accounts for the magnetic-field vacuum fluctuations which cannot be neglected in the case 
of metallic particles. Finally, the effect of finite temperature can be easily addressed. We 
note that a different approach has been recently presented^ where the EM Green's ten- 
sor entering the fluctuation-dissipation theorem is calculated by means of a finite-difference 
frequency-domain method. 

When a particle is a member of a colloidal crystal and a net vdW force exerted on the 
particle is evident (e.g., in a finite slab of a colloidal crystal or in an infinite crystal containing 
point and/ or line defects), it is calculated from a pairwise addition of the forces stemming 
from the all the other particles of the crystal. So, at first glance, an extension of the Ref. [17 
to the case of a colloidal system would be based on a pairwise summation of the (exact) 
force for a pair of particles. However, such an approach is only approximatel y co rrect since 
the vdW interactions are not additive. The way to extend the method of RefJl7l to the case 
of a colloidal crystal is to derive a semi-analytical expression of the EM Green's tensor for 
the particular crystal. The knowledge of the EM Green's tensor everywhere in space allows 
the calculation of the cross-spectral correlation functions of the vacuum EM field which are 
contained in the EM Maxwell stress tensor, by application of the fluctuation-dissipation 
theorem. By integrating the Maxwell stress tensor over the surface of the particle we obtain 
the vdW force. The paper is organized as follows. In section [IT] we provide a brief overview 
of fluctuational electrodynamics and the Maxwell stress tensor. In section II II I we provide 
expressions for the EM Green's tensor, (a) for arbitrary collections of a finite number of 
scatterers, (b) for infinite, periodic and defected crystals, and (c) for finite slabs of colloidal 
crystals. In section [IV] we apply the formalism to the case of a monolayer of polystyrene 
spheres containing a single defect. Section [V] concludes the paper. 

II. VAN DER WAALS FORCE 
A. Maxwell stress tensor 

We consider a finite scatterer with electric permittivity e s and/or magnetic permeability 
fA s different from those, e^, fih of the surrounding homogeneous medium. According to 
classical electrodynamics, the exerted force F on a finite scatterer in the presence of electric 
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E and magnetic H fields satisfying the Maxwell equations is obtained by integrating the 
time-average Maxwell stress tensor T,y^ over the surface around the scatterer 

(F t ) t = [ J2(^MidS (1) 

where (..)t denotes the time average, n is the normal vector at the surface surrounding the 
object, and i,j = x,y,z. The components of the tensor (Tij) t are given by 

(Ty)i = e h e (Ei(r,t)Ej(r,t))t + fi h fi (H i (r,t)H j (r,t)) t 

~ 2 5 v ^ ( r ' ^ Ei ' ( r ' *)>* + VhVo ( Hi> ( r ' *) H * ( r ' *)>*! " ( 2 ) 

v v 

Sij is the Kronecker symbol and eo, /io are the electric permittivity and magnetic permittivity 
of vacuum, respectively. 



B. Fluctuation-dissipation theorem 



In the absence of other radiation sources, the fields E, H are generated by the thermal 
radiation emitted from the same or neighboring scatterers at finite temperature (thermal 
fluctuations) or by vacuum radiation at zero temperature (zero-point fluctuations). The 
time-correlation function (Ei(r,t + r)Ej(r' ,t))t contained in Eq. ([2]) is calculated within the 
framework of fluctuational electrodynamics,- 1 ^ namely iromM 



(E i (v,t + r)E j (r',t)) t = Re 



J 2^ exp ( luJT > W > 



EE i 
ij \ 



(3) 



The quantity (r, r'; to) is the cross-spectral correlation function for the electric field. 
For a system at thermal equilibrium, i.e., the scatterer, the surrounding medium and its 
neighbouring scatterers at the same temperature T, Wij is provided by the fluctuation- 
dissipation theorem^ 1 ^ 

1 



W? E {v,r';u) 



4a;/i/ t /ioC 2 ImG^ £ '(r, r'; uj) hio 



1 + 



(4) 



where h is the reduced Planck's constant, kg is the Boltzmann's constant and Gfj (r, r'; to) 
is the component of the full Green's tensor Gij which provides the electric field at r due to 
an electric dipole source at r'. The time-correlation function (Hi(r,t + r)Hj(r f , t)) t for the 
magnetic field is given similar to Eq. ([3]) with W EE substituted by 

1 



W? H (r, r'; u) = Auue h e (?\mG* H {r, r'; u)tu 



1 + 



exp(hu/kBT) — 1 



(5) 



We note that, the final value of the vdW force acting on a scatterer is obtained by subtracting 
from Eq. ([1]) the force which remains in the absence of the scatterer as it is the case for the 
calculation of the Casimir force between two semi-infinite slabs.— However, in vacuum, the 
Green's tensor and the corresponding Maxwell stress tensor, Eq. (j2j), are constant in space 
and their integral over a closed surface is zero. From the above, it is obvious that the central 
quantity which essentially determines the force acting on the scatterer is the EM Green's 
tensor. 

III. ELECTROMAGNETIC GREEN'S TENSOR 

A. Multipole expansion of the EM field 

Let us consider a harmonic EM wave, of angular frequency to which is described by its 
electric-field component 

E(r, t) = Re [E(r)exp(-iwt)] . (6) 

In a homogeneous medium characterized by a dielectric function e(uj)eo and a magnetic 
permeability fi(u)fio, where eo, fio are the electric permittivity and magnetic permeability of 
vacuum, Maxwell equations imply that E(r) satisfies a vector Helmholtz equation, subject 
to the condition V • E = 0, with a wave number q = u/c, where c = 1/^/ze/ioeo = co/v^ 
is the velocity of light in the medium. The spherical-wave expansion of E(r) is given by^ 

E M = J2 \ a Eimfi(qr)X lm (r) + a Elm -V x [/,(gr)X, m (r)] I , (7) 

1=1 m =-i ^ Q J 

where a[ m (P = E,H) are coefficients to be determined. X/ m (f) are the so-called vector 
spherical harmonics^ and fi may be any linear combination of the spherical Bessel function, 
ji, and the spherical Hankel function, . The corresponding magnetic induction, B(r), can 
be readily obtained from E(r, t) using Maxwell's equations^. 

B. Scattering from a single scatterer and the corresponding Green's tensor 

In this subsection we present a brief summary of the solution to the problem of EM 
scattering from a single sphere (Mie scattering theoryi^^) along with the expression for the 
single-sphere Green's tensor. We consider a sphere of radius S, with its center at the origin 
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of coordinates, and assume that its electric permittivity e s and/or magnetic permeability 
fi s are different from those, eh, \Xh of the surrounding homogeneous medium. An EM plane 
wave incident on this scatterer is described, respectively, by Eq. ([7|) with // = j t (since the 
plane wave is finite everywhere) and appropriate coefficients a° L , where L denotes collectively 
the indices Plm. That is, 

E°(r) = $>°J L (r) (8) 

L 

where 

Jeim(r) = —V x ji(q h r)X lm (i), J H im{r) = ji(q h r)Xi m (r) (9) 

Qh 

and qh = t^ii^ui j c§. The coefficients a° L depend on the amplitude, polarization and propa- 
gation direction of the incident EM plane wave and are given by Eqs. ( 1371) (subsection IIII Gil 
for g = 0. 

Similarly, the wave that is scattered from the sphere is described by Eq. (J7J) with 
fi = , which has the asymptotic form appropriate to an outgoing spherical wave: 
hf ~ (— i)' exp(i<27ir)/ig/jr as r — > oo, and appropriate expansion coefficients a\. Namely, 

E + (r) = 5>+H L (r) (10) 

L 

where 

H Eim (r) = —V x hf(q h r)X lrn (r), H Hlm (r) = hf (q h r)X lm (r) . (11) 

Qh 

The wavefield for r > S is the sum of the incident and scattered waves, i.e., E out = E° + 
E + . By applying the requirement that the tangential components of E and H be continuous 
at the surface of the scatterer, we obtain a relation between the expansion coefficients of the 
incident and the scattered field, as follows: 

V 

where Tiy are the elements of the so-called scattering transition T-matrix.— Eq. (fl2j) is 
valid for any shape of scatterer; explicit relations of the T-matrix for scatterers of various 
shapes can be found elsewhere.— 1 ^ 

The Green's tensor for a single sphere is given by 25 

G«(r,rO = -iu,^!^ (13) 
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The vector functions Ri-i(r), Rm(r) are dimensionless eigenfunctions of the wave operator 

A(r) = -j^prV x V x (14) 
e(r)/i(r) 

for a single scatterer which are regular at its center.— ^ The vector functions Il-^t), 
are also eigenfunctions of the operator ( JT4l) but they are infinite at the sphere center.— ^ 
The Green's tensor of Eq. (TT3l) will be the basis for the construction of the corresponding 
tensor for a collection of spheres. 

C. Green's tensor for many scatterers 

We consider a collection of N nonoverlapping scatterers described by a permittivity e s 
and permeability fi s centred at sites R n in a homogeneous host medium described by e^, fih, 
respectively. In site-centered representation, the Green's tensor for the system of scatterers 
satisfies^ 6 . 

^ [uj 2 5i"i - Aj»i(Rn + r n )] Gii'(R„ + r n , R n / + r' n /) = <W<Kr n - r' n /)<W (15) 

i 

where r n = r — R n , r' n / = r' — R n /, and i, %' = x,y,z. The operator Aj«j(r) is given by 
Eq. dHJ). It can be verified that the Green's tensor satisfying Eq. (IToT) is the following—^ 

G ii ,(R n + r n ,R n , + T' n ') = G$ n (r n ,r' n ,)5 nn/ -i J 6h ^ 3 * ^Ki( r n) D uLR n u-,A r 'n')- ( 16 ) 

° LU 

G^) n (r n ,r' n /) is the Green's tensor for a single scatterer located at R ra and it is given by 
Eq. (ITB1) . The vector functions i?£.j(r n ), Rui( r n) are the dimensionless eigenfunctions of the 
operator of Eq. (|T4|) for the sphere at R n . -D22' are propagator functions that represent the 
contributions of all possible paths by which a wave outgoing from the n'-th scatterer produces 
an incident wave on the n-th scatterer, after scattering in all possible ways (sequences) by 
the scatterers at all sites including the n-th and n'-th scatterers. The specific form of the 
-D^z/ propagator functions depends on the geometrical arrangement of the scatterers. 

D. Propagator for an arbitrary collection of scatterers 

For an arbitrary collection of a finite number N of scatteres, the D-propagator is given 
by^ 

nmi' r\nn' , \ ^ \ ^ \ ^ T~\nn" rpn" r)n"n' ( -I ij\ 

U LU — lL LU ' 2-^i / j LL" L"U" 1 L L"'U ' \ LI ) 

n" L" U" 
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The matrix f^L' appearing in Eq. (ITT)) is called free-space propagator and transforms an 
outgoing vector spherical wave about R n / in a series of incoming vector spherical waves 
around R n -— The matrix T£ L , is the scattering T-matrix of a scatterer of general shape,— ^ 
located at R n . 



E. Propagator for periodic arrays of scatterers 

For the case of an infinite number of same spheres arranged periodically, in one- (ID), 
two- (2D) or three (3D) dimensions, the propagator D^jj * s gi ven as a Fourier transform 

D% = - [ d^expfik ■ (R n - R n ,)]D LL ,(k), (18) 
v Jbz 

where q is the space dimensionality, the integration in Eq. ( fl8l) is carried out within the 
Brillouin Zone (BZ), k is the Bloch wavevector, and v is the BZ volume. R n are the Bravais 
lattice vectors. D LL /(k) is given by 

D LL ,(k) = n LL ,(k) + D LL ,,(k)T L//L/ M L/ ,, L/ (k). (19) 

L"L"> 

Tl"l"' is the T-matrix of the spheres. fiz,L'(k) depend only on the crystal lattice and 
are known as structure constants a term which is common in the Korringa-Kohn-Rostoker 
method^ for the calculation of the electronic band structure of atomic solids. They can be 
found by Ewald-summation techniques.— 1 ^ Eqs. (Hj) and (jSJ) require the calculation of the 
Green's tensor [via Eq. ( IT5I) ] for an infinitely periodic lattice of scatterers; therefore, only 
the D L ° L , component (that for the central unit cell) is needed since all spheres are equivalent 
for the case of a Bravais lattice with one sphere per unit cell. 

We note that the propagator of Eq. ffT9]) does not yield a net nonzero vdW force, since 
it corresponds to an infinitely periodic system. However, the propagator of Eq. (fl9l) can be 
used as a basis for calculating the corresponding propagator of a system containing, e.g., 
one or more point defects (not symmetrically distributed within the crystal), in which case 
a net vdW force emerges. If, for example, the colloidal particles (described by a a scattering 
matrix Tq LU ) positioned at R n in an otherwise periodic crystal, are substituted by other, 
different particles, each of them described by a scattering matrix T£ L ,, the propagator of the 
defected system is given similar to Eq. ( ITTj) . i.e., 

Dit = D n 2' L , + E E E Dii;,ATz;; L ,„D&; L , (20) 

n" L" V" 
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where AT^" L „, = T£,'' L „, — T^ LL , and -Dqll' * s ^ ne propagator of the periodic system given by 
Eqs. (fUD and (fT9|). 

F. Propagator for finite slabs 

In reality, the colloidal systems are not infinitely periodic but they are actually slabs 
consisting of a finite number of planes of particles (scatterers). In this case, the vdW force 
exerted on a given scatterer depends on the position of the plane within which it is located 
and can therefore be very different for a scatterer on a surface plane than a scatterer at an 
innermost plane. In the following lines, we will provide a formalism for the propagator for a 
slab consisting of N p planes of scatterers. It is assumed that all the planes of the slab have 
the same 2D periodicity with the associated lattice vectors given by 

R n = niai + n 2 a 2 , (21) 

where ai and a 2 are primitive vectors in the xy plane and ni, n 2 = 0, ±1, ±2, ±3, ■ ■ -. The 
corresponding 2D reciprocal lattice is defined by 

g = mibi + m 2 b 2 (22) 

where mi, m 2 = 0, ±1, ±2, ±3, • • ■ and bi, b 2 are primitive vectors defined by 

hi ■ aj = 2ir5ij , i,j = 1, 2. (23) 

Although each plane of the slab must have the same 2D periodicity, the spheres within each 
of the N p planes can be different in terms of shape, size or refractive index. 

The propagator for a scatterer residing at the v-th plane [y = 1, 2, • • • , N p ) of a slab is 
written as a sum of three terms^ 

F u^LU = D u%L' + ^2 ^2 P ^LL» T v;L"L»>D^L"'L> + P u^LL' ( 2 4) 
n L" L"> 

The matrix D™™ L , represents all the possible scattering paths within the v-th plane by 
which a wave outgoing from the m-th sphere of this plane produces an incident wave on 
the n-th sphere of the same plane, after scattering in all possible ways by all the spheres of 
this plane including the central sphere (every sphere represented by the scattering matrix 
T v -lv)- It is given by application of Eq. ffl8|) to a 2D periodic lattice, i.e., 

9 



DTlV = J" / / eX P( ik ll ■ Rn m )A,;LL'(k||) (25) 

>jo 7 J SBZ 

where 

D^(k|,) = ^[[i-nckii)^]- 1 ] n L » L ,(kn) (26) 

L" 

where R nm = R n — R m , So is the area of the Surface Brillouin Zone (SBZ) corresponding 
to Eq. (1221) . and f2^/(k||) are the 2D structure constants. 

The matrix P® n LL / appearing in the second and third terms of Eq. ( l24|) represents all 
scattering paths by which an outgoing wave from the n-th sphere of the v-th plane exits 
from that plane to produce an incident wave on the central sphere of the same plane after 
scattering in all possible ways by all the planes of spheres of the slab, including the u-th 



plane. In the next subsection we wi 
F®? LL , which is given in detail in Ref. 



lpresent a summary of the derivation of P^^i and 
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G. Calculation of P°™ L , and F™ LL , 

A wave outgoing from the n-th sphere of the v-th plane has the form of Eq. (flOj) 

E"(r) = ^6+(n;i/)H L (r) (27) 

L 

where r nu is the position vector with respect to the center of the n-th sphere of the u-th 
plane. We can expand the wave of Eq. (1271) into a sum of plane waves propagating or 
decaying away from the u-th plane as follows.— To the right of the u-th plane we have 



E aat +/ x 



IT [ f ^llS E r' + (k||)exp[iK+-(r-A 2 (z/))] (28) 

bo J J SBZ „ 



with 

Eg + (k,|) = exp[-i(k|, R n - K+ d 2 (i/))] J2 A L;i (K+)6+(n; „) (29) 

L 

where i = 1,2. A 2 (z/) is a reference point on the right of the u-th plane at d 2 (u) from its 
center (see Fig. [I]). To the left of the ^-th plane we have 



B out - (r) 



x - f f d% { K ut "( k n) ex p[ iK § • ( r - M"))) ( 3 °) 

3 J J SBZ „ 
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with 

E% "(k||) = exp[-i(k|| R n + K g d x {v))\ A i;i (K-)6+(n; i/) (31) 

L 

where Ai(z/) is a reference point to the left of the u-ih plane at — d\(v) from its center 
(see Fig. [1]). Kg is given by Kg = I ky + g, ± q 2 — (ky + g) ) , where the +, — sign 
defines the sign of the z component of the wavevector. The coefficients A^.j are given from 
Eqs. (19) and (20) of Ref . \31 



The plane waves of Eq. ( 128]) will be multiply reflected between two parts of the slab, the 
first (right part) consisting of all planes to the right of the u-th plane, and the second (left 
part) consisting of all planes to the left of the (v + l)-th plane (including the u-th plane), to 
produce a set of plane waves incident on the v-th plane from the right, which we can write 
formally as follows 



«~>i 



l - f f d 2 fc||^E--(k||)exp[iKg .(v-A 2 (u))] (32) 

3 J JSBZ „ 



with 



e% "(k B ) = £ (Q m (^ 2 )[ J - Q> + !; 1 )Q III ( z/ ' 2 )1 _1 W ^' + ( k ll) ( 33 ) 



where Q 1) and Q IH (z/; 2) are the appropriate matrices which determine the reflection 

(diffraction) of a plane wave by the left and the right parts of the slab respectively, as defined 
above. These matrices are shown schematically in Fig. [TJ 

Similarly, the plane waves of Eq. (I3"U|) will be multiply reflected between two parts of the 
slab, the first (left part) consisting of all planes to the left of the v-th plane and the second 
(right part) consisting of all planes to the right of the [y — l)-th plane (including the ^-th 
plane), to produce a set of plane waves incident on the v-tk\ plane from the left, which we 
can write formally as follows 



E' 



1 / / d\\ E ? + ( k ll) ex Pt iK ^ ■ ( r " A iW)] (3 4 ) 

^0 J JSBZ „ 



with 



*S + ( k n) = E (Q n (^ - Q 111 ^ - 1; 2)Q n ("; i)]" 1 }^ ^"( k n) (35) 



where Q n (z^; 1) and Q HI (z/ — 1; 2) are again the appropriate matrices, shown schematically 
in Fig. [U A more detailed description of these matrices and the way these are calculated is 
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to be found in Ref. [3l|. We note that for v = 1(N) we have only waves incident from the 
right (left). 

Each plane wave in Eqs. (1321) and (1341) can be expanded in spherical waves about the 
central sphere of the v-th plane in the manner of Eqs. ([H]) and For a plane wave 

coefficients are given by^ 



v k||)exp[iK g • (r — A 2 (V))], incident on the z/-th plane from the right, the multipole 



a°(K g ) = exp[-iK g • d 2 (u)) £ A° ;4 (K g )E* "(k,. 



(36) 



And for a plane wave, E g ™ + (k||) exp[iK g • (r — Ai(i/))], incident on the z/-th plane from the 
left, the multipole coefficients are^i 



a° L (K+) = exp[iK+ ■ d^)] £ A%^)E™ +(k„ 



(37) 



where AP LH are given by Eqs. (12) and (13) of Ref. [3l|. 

Finally, to obtain the wave incident on the central sphere of the z/-th plane, which derives 
from the outgoing wave of Eq. (1271) . we must add to the waves given by Eqs. (132]) and (|34|) 
that which is due to the wave scattered from all the other spheres of the v-t\i plane and it 
is given by multiplying the coefficients a° L of Eqs. ( I36p and (l3"Tj) by the multiple-scattering 

for the z/-th plane of spheres. We have 



matrix ^[1 - QT^ 1 
v 



LL' 



1 

$0 J JSBZ 
1 




g s=± V 



J LL 



Sn 




d 2 k\\ exp(— iky • R n ) 



SBZ 



(38) 



[I-ftTJ- 1 ^ 



LL 



(39) 



where T u: ll' is a matrix defined by 

r,.p JwM w(k||;w) = ^^|exp[-i(K g - K+) • d 2 (v)]^ m;i (K 

g>* g',j' 

x [Q m (z,;2)[I-Q II (z, + l;l)Q III (z,;2)]- 1 
+ exp[i(K+ - K g ,) • d 1 (v)]A^ ni . i (K+) 

x Q I V;i)[i-Q m (^-i;2)Q I >;i)]- 1 



Ap'I' m ' ; i'(K . 
gi;g'i' & 



(40) 
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Therefore, from Eq. (|39|) . P„ n LL i is given by 

P» n LL> = 4- / / rf2fc ll exp(-ik|| ■ R n ) [[I - nT v ]- l T v ] (41) 

•JO J J SBZ 1 J Li ' 

Accordingly the second term in Eq. §M§ becomes^ 

E E E P^ L v>T v ,l^D:1, iil , = i- / / _ d 2 ^ [[I - nTj-^T^] ^ (42) 



„ L „ L ,„ S J J SBZ 



where D u . LL f(\i\\) is given by Eq. fl26|) . Finally, the matrix F®° LL ,, defined by Eq. (|24"|) . is 
given by 



[i - ot,]- 1 ^ + r„(i + t v [i - fiT^-'n)} (43) 

LL 1 



F ™LL> = f / / rf^ll 

IV. NUMERICAL EXAMPLE 

The evaluation of the propagator, either from Eq. (125]) or Eq. (|43|) . requires a numerical 
integration over the entire SBZ. Using symmetry to reduce the area of integration to a part 
of SBZ is not profitable in the present case. However, when one deals with scatterers whose 
dielectric function contains a positive imaginary part, the intergrand in Eqs. (121)1) or fT4"3"|) 
is a relatively smooth function of ky , and the integration can be performed without much 
difficulty by subdividing the SBZ (a square in our example) into small squares, within which 
a nine-point integration formula^ is very efficient. Using this formula we managed good 
convergence with a total of 576 points in the SBZ. 

When computing the vdW for T = 0, we first integrate the Maxwell stress tensor for 
a specific frequency over the surface of the body and afterwards we perform the frequency 
integration, i.e., the vdW force F is calculated by integrating the force spectrum F(u): 
F = J °° F{uj). Both integrals are obtained numerically. We note that, in the Lifshitz theory 
for half-spaces,— the frequency integration is done analytically using contour integration. 
The numerical integral over frequencies is convergent since, in the limit of uj — > oo, the 
refractive index of most materials tends to unity and the corresponding Green's tensor of 
the system tends to that of vacuum which is constant in space. However, the integral over 
a closed surface of a constant tensor vanishes and therefore F(uj) — > as uj — > oo. 

We consider the case of a 2D square lattice (monolayer) of polystyrene nanospheres of 
radius 10 nm. The dielectric function of the spheres which is generally complex for high 
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frequencies is taken from numerical fit to experimental data.- We have calculated the force 
acting on a single polystyrene nanosphere when we remove one of its first neighbouring 
spheres (see inset of Fig. [2]). In this case, we first calculate the propagator for the periodic 
square lattice from Eqs. (|25p and f[2Ej) (which yields vanishing net vdW force) and then make 
use of Eq. (T2Uj). 

In Fig.[2]we show the net vdW force (x- and y- components) for different lattice constants 
a of the the underlying 2D square lattice. While each of the components oscillates from 
positive to negative values, it is evident that there exists a value of the lattice constant a, 
namely a ~ 47 nm, where the net force is zero and this particular sphere rests in equilibrium. 
Overall, the magnitude of the vdW force decreases with the lattice constant, as expected. 

V. CONCLUSION 

We have presented a method for the calculation of the vdW forces in colloidal systems 
such clusters of colloidal particles, infinite periodic or defected crystals, and colloidal crystals 
slabs. The method is based on the fluctuation-dissipation theorem which relates the cross- 
spectral correlation functions entering the formula for the vdW force (integral of the Maxwell 
stress tensor over the particle surface) with the EM Green's tensor of the system of particles 
(scatterers). The calculation of the Green's tensor is based on a rigorous multiple-scattering 
formalism for EM waves. The accuracy stems from the fact that it does not include any kind 
of approximations apart from the unavoidable cutoffs in the angular momentum expansion 
and/ or in the plane- wave expansion of the EM field. As such, the method includes all 
essential multipole terms beyond the dipole term in the EM response of the scatterers and 
is valid for any distance between the scatterers. By including a priori all the possible 
multiple-scattering processes of the vacuum fluctuations, the method, naturally, accounts 
for all possible many-body interactions between the scatterers and therefore goes beyond 
the approximation of pairwise interactions. 

Finally, we note that a theoretical approach, analogous to the multiple-scattering treat- 
ment for the wave equation, has been developed for solving the Poisson equation in solids 
described by arbitrarily shaped, space-filling charges.— By combining this electrostatic 
multiple-scattering approach with the vdW theory presented in this work, one can devise a 
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general, first-principles theory for the determination of colloidal structure. 
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FIG. 1: The Q-matrices appearing in Eq. (|40|) . The position vectors di, d2 of the z^-th layer along 
with the corresponding origins Ai, A2 are also shown. 



16 



-3.5 1 1 1 1 1 1 1 1 1 1 

20 30 40 50 60 

Lattice constant (nm) 

FIG. 2: (Color online) Inset: 2D square lattice of 10 nm polystyrene spheres containing a single 
defect (one missing sphere). Graph: the x- (squares) and y- (circles) component of the vdW force 
exerted on a single polystyrene nanosphere when its right neighboring sphere is missing. 
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